Universal amplitudes in the FSS of three-dimensional spin models 
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In a MC study using a cluster update algorithm we inves- 
tigate the finite-size scaling (FSS) of the correlation lengths of 
several representatives of the class of three-dimensional classi- 
cal 0(n) symmetric spin models on the geometry T 2 x R. For 
all considered models we find strong evidence for a linear re- 
lation between FSS amplitudes and scaling dimensions when 
applying antiperiodic instead of periodic boundary conditions 
across the torus. The considered type of scaling relation can 
be proven analytically for systems on two-dimensional strips 
with periodic be using conformal field theory. 

PACS numbers: 64.60.Fr, 75.10.Mk, 75.40.Mg, 11.25.Hf 



Conformal invariance of 2D systems at a critical point 
has turned out to be the key feature for a complete, ana- 
lytical description of their critical behavior . In par- 
ticular, conformal field theory (CFT) supplies exact FSS 
relations including the amplitudes for these 2D models. 
For strips of width L with periodic boundary conditions, 
i.e. the S 1 x R geometry, Cardy has shown that the 
FSS amplitudes of the correlation lengths £j of primary 
(conformally covariant) operators are entirely determined 
by the corresponding scaling dimensions Xf. 
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with a model independent overall amplitude A = 1/2tt. 
This result relies on both, the greater restrictive strength 
of the 2D conformal group compared with the higher di- 
mensional cases, which is needed for the definition of the 
"primarity" of operators, and the fact that the considered 
geometry is conformally related to the corresponding flat 
space R 2 . 

Generalizing these results to more realistic 3D geome- 
tries within the CFT framework generically destroys the 
rich 2D group structure. Keeping at least the conformal 
flatness condition, Cardy Q arrived at a conjecture of the 
form (|l]) for the S' n_1 x R, n > 2 geometries. Mainly for 
reasons of the numerical inaccessibility of these geome- 
tries Hcnkcl H^j considered the situation where even this 
latter condition is cancelled: investigating the scaling be- 
havior of the S = | Ising model on 3D columns T 2 xl 
with periodic (pbc) or antiperiodic (apbc) boundary con- 
ditions across the torus via a transfer matrix calculation, 
he found for the correlation lengths of the magnetization 
and energy densities (the only primary operators in the 
2D model) in the scaling regime the ratios: 



£<r/£ e = 3.62(7) periodic be, 
^a/S,e = 2.76(4) antiperiodic be. 



(2) 



Comparing this to the ratio of scaling dimensions of 
x t /x a — 2.7326(16) a relation of the form (|l|) seems not 
to hold, unless the boundary conditions are changed to be 
antiperiodic. This is in qualitative agreement with nu- 
merical work done by Weston H. 

In this letter, we first revisit the Ising model on the 
T 2 x R geometry trying to decide the exposed question 
with an independent Monte Carlo (MC) method and at 
an increased level of accuracy. The main purpose is to 
investigate further models - in our case 0(n), n > 1 spin 
models -, thus adding evidence that Henkcl's result is 
not just a numerical "accident" but reflects a universal 
property of such 3D systems. 

The model — We consider an 0(n) symmetric classi- 
cal spin model with nearest-neighbor, ferromagnetic in- 
teractions in zero field with Hamiltonian 



n = -jJ2 



3 j< 
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The spins are located on a sc lattice of dimensions 
(L x , L y , L z ) with L x — L y , modeling the T 2 geometry 
by applying periodic or antiperiodic be along the x- and 
y-directions. Effects of the finite length of the lattice in 
the z-direction are minimized by choosing L z such that 
Lz/S, ^ 1 and sticking the ends together via periodic be. 
As is well known all of these models undergo a con- 
tinuous phase transition in three dimensions, so that at 
the critical point the correlation length diverges linearly 
with the finite length L = L x . Particular representatives 
of this class are the Ising (n = 1), the XY (n = 2) and 
the Heisenberg (n — 3) model. 

The simulation - For our MC simulations we used 
the Wolff single-cluster update algorithm [0] which is 
known to be more effective than the Swendsen-Wang |ic[ | 
update for three-dimensional systems As we want 

to consider antiperiodic be for all systems in addition 
to the generic periodic be case, the algorithm had to be 
adapted to this situation using the fact that in the case of 
nearest-neighbor interactions antiperiodic be are equiva- 
lent to the insertion of a seam of antiferromagnetic bonds 
along the relevant boundary. 

The primary observables to measure are the connected 
correlation functions of the spin and the energy density: 



G°(xi,x 2 ) 
G?(xi,x a ) 



(s(xi) • s(x 2 )) 
(e( Xl )e(x 2 ))- 



-<s)<s> 
(*><e>. 



(4) 
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The correlation lengths & in Eq. ([l]) being understood 
as measuring the correlations in the longitudinal In- 
direction, one may average over estimates G c (xi,X2) 
such that (xi — X2) || e z and i = |xi — X2I = const, 
thus ending up at estimates G c '"(i). This average can be 
improved by considering a zero momentum mode projec- 
tion |l2] ], i.e., by correlating layer variables made up out 
of the sum of variables in a given layer z = const instead 
of the original spins or local energies; this reduces the 
variance by a factor of 1/L 2 ., the influence of transversal 
correlations being irrelevant for large distances i p3| . 

Assuming an exponential long-distance behavior of 
the correlation functions @), extracting the correlation 
lengths via a straightforward fitting procedure requires a 
nonlinear three-parameter fit of the form 



qc,\\ (j) = gcII(o) exp (-»/£) + const, 



(5) 



since any numerical estimation of G c, "(i) necessarily fails 
to reproduce the correct long distance limit G c, "(i) — > 
as i — > 00 exactly. As this amounts to an investment 
of the gathered statistics into the determination of three 
parameters, two of which are completely irrelevant for 
our ends, we used an alternative method which intrinsi- 
cally eliminates the two irrelevant parameters by using 
differences and ratios of G c '"(i) rather than the values 
themselves. Given the correlation function behaves as 
(pf) , estimators £j for the correlation length are given by: 
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G c 'W(i) - G c ^(i - A) 
G c 'H(i + A) - G c 'H(i) 



(6) 



The generic value for A is one, but it might be advan- 
tageous to choose A > 1 in order to enhance the local 
drop of G c '"(i) between i and i + A (the signal) against 
the fluctuations (the noise) . Following this procedure one 
ends up with a set of estimators for the correlation length 
as a function of distance i as depicted in Fig. [I] for the 
spin-spin correlations of the Ising model: after a transi- 
tion regime starting at i = A which is a consequence of 
the discreteness of the lattice as well as the above men- 
tioned zero momentum mode projection, the estimates 
settle at a plateau indicating that the exponential long 
distance behavior has been reached. 

The error bars in Fig. [l] were generated using a com- 
bined binning and "jackknife" resampling scheme [ |l4||T5| 
which is necessary due to the strong non-linearity of the 
transformation (|6|); on the same grounds we checked for 
the necessity of a bias correction. Final values for the 
correlation lengths were obtained by an average of the 
estimators £j ; as neither the estimates for very small dis- 
tances i nor - because of the periodicity of the lattice in 
the z-direction - those for distances i > L z /2 are reliable 
estimates for the continuum correlation length, the range 
of distances i m in, ■ ■ • , 'max to average over was determined 
by a procedure of statistical optimization, a generalized 



% 2 -test. In order to minimize the theoretical variance of 
the final average £ each element £j was weighted by a 
factor proportional to a row sum of the inverse covari- 
ance matrix, this matrix itself being again estimated by 
a jackknife technique ]13|,^6[. 




FIG. 1. Example of the set of estimators £i for the mag- 
netization density of a 18 2 x 214 Ising system with periodic 
be. The typical distance A in Eq. (^|) was set to 8. 

The Ising model - Simulations of the Ising model 
were done at the most accurate estimate for the bulk in- 
verse critical temperature available, (3 C — 0.2216544(3) 
Jl7t , where the influence of the given error in (3 C on the 
results for the correlation lengths was checked via a tem- 
perature reweighting technique and found negligible com- 
pared to the statistical errors; this applies to the other 
models considered in this note as well. To be able to per- 
form a FSS analysis, simulations were done for system 
sizes between 4 2 x 48 and 30 2 x 356 rs 3 x 10 5 sites, accu- 
mulating about eight million independent measurements 
for each system. 

As is obvious from the example in Fig. ^(a) the final 
estimates for the correlation lengths show up an almost 
perfect linear scaling behavior as a function of the trans- 
verse system size L x . Considering the amplitudes £,/L x 
reveals, however, that corrections to the leading linear 
scaling behavior are relevant and can be clearly resolved 
within the accuracy of the data, cp. Fig. ||(b). In order 
to extract the leading amplitudes in the scaling regime 
nonlinear fits of the form 



£(L X ) = AL X + BL\ 



(7) 



were done. Even though some field theoretical estimates 
for the correction exponents exist || , we decided to keep 
a as a parameter, ending up at an effective correction ex- 
ponent that takes higher order corrections into account, 
which have some importance for the small systems; suc- 
cessively dropping systems from the small L x end while 
monitoring the goodness of fit parameters \ 2 and Q then 
acts as a consistency check. As a rule, the overall cor- 
rections are negative for systems with periodic be and 
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positive in the case of antiperiodic be. 

As a result of this fitting procedure we arrive at the 
following final estimates for the amplitudes A in Eq. (]?]) 
and their ratios: 



A a = 0.8183(32) 
A e = 0.2232(16) 
A a j A e = 3.666(30) 

Aa = 0.23694(80) 
A e = 0.08661(31) 
A a J A e = 2.736(13) 



for periodic be, 



for antiperiodic be. 



(8) 



Comparing this to the ratio of scaling dimensions 

BIB 



(1 - a)/v _ 2{vd - 1) 
vd 



= 2.7326(16), (9) 



we find that the amplitude and exponent ratios agree very 
precisely in the case of antiperiodic be across the torus, 
while in the periodic case they differ by an amount of 
some thirty sigma. In comparison to a first exploration 
by Weston 0, who found ratios of about 3.7 for peri- 
odic and 2.6 for antiperiodic be, the precision could be 
increased by over an order of magnitude. 
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FIG. 2. (a) FSS plot for the spin correlation length £cr(L x ) 
of the 3D Ising model with antiperiodic be. (b) Scaling of 
the amplitudes £<jJ_Lx- Solid lines represent least-square fits 
according to Eq. 



XY and Heisenberg model - Although being strin- 
gent in itself, up to this point the above result is a sin- 
gular, maybe casual, statement for the special case of 
the Ising model. Believing in a universal law needs a 
broader backing with successful examples, two of which 
being considered here. 

Simulations for the XY and Heisenberg model were 
done at the estimated inverse critical temperature values 
(3 C = 0.4541670(32) and (3 C = 0.693004(7), respectively, 
which are weighted means of recent literature estimates 
pd] , ^9| -[2l]] . Investing about three years of workstation 
time for both models altogether and using the same sys- 
tem sizes as in the Ising case, we took between four and 
eighteen million independent measurements for each sys- 
tem, both for periodic and antiperiodic be. Applying the 
outlined tools of data analysis we arrive at scaling and 
amplitude plots similar to those in Fig. ^[ Traversing 
the above described fitting procedure leads to final es- 
timates for the amplitudes A a and A e according to Eq. 
([?]), which are shown in Table |. Comparing the results 
for the ratios A a /A € with the ratio x t jx a of scaling di- 
mensions, we arrive at a highly precise agreement for the 
case of antiperiodic be and an obvious divergence in the 
standard periodic be situation for both, the XY and the 
Heisenberg model. Thus a linear relation between scal- 
ing amplitudes and scaling dimensions according to Eq. 
([j]) is almost certainly valid for three generic, non-trivial 
examples of 3D spin models, and one might well assume, 
that it is satisfied for the whole class of O(n) spin models, 
a view which is supported by further simulations for the 
n = 10 case [|16[ and an analytic result for the limiting 
case n — > oo |22j ], which is known to be equivalent to the 
spherical model |2^] . In view of the analogous 2D results 
it is not too far fetched, then, to argue that the numerical 
results provide evidence that this relation might be of a 
universal, model independent kind. 
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FIG. 3. Amplitudes A according to Eq. (|ic|) versus the in- 
verse dimension of the order parameter 1/n for n = 1, 2, 3, oo. 

Universal amplitudes - Given the scaling ampli- 
tudes for the 3D systems with antiperiodic be behave 
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according to Eq. (pi), one may ask further, what the am- 
plitude A in Eq. (]!]), that was l/2n in the 2D case, be- 
comes in the 3D scenario and, furthermore, if it holds 
true that it is universal in the sense that all the model- 
dependent information is condensed in the scaling dimen- 
sions Xi. The transfer matrix approach cannot give an 
answer to this question, because in the Hamiltonian limit 
amplitudes are only given up to an overall normalization 
factor. Supposed such a relation holds, from the above 
results one can give an estimate for these amplitudes us- 
ing the amplitudes A a , which are usually more accurate 
than A e . Using the scaling dimensions of the spin of 
x a = 0.5175(5), x a = 0.5178(15) and x a = 0.5161(17) 
for the Ising, the XY and the Heisenberg model, respec- 
tively, one has: 

f 0.12262(43) Ising 
A = A a x a = < 0.12486(47) XY . (10) 

[ 0.12625(49) Heisenberg 

Taking into account the corresponding amplitude of 
the spherical model, which is A « 0.13624 l^p^l, and 
comparing the variation of these values with the given 
errors, as is shown in Fig. ||, it becomes clear that these 
amplitudes do in fact depend on the model under con- 
sideration and seem to vary smoothly and monotonically 
with the dimension n of the order parameter. 

To summarize, the amplitudes of the FSS of the corre- 
lation lengths of the magnetization and energy densities 
of O(n) spin models are linearly related to the corre- 
sponding scaling dimensions for the T 2 x R geometry 
when choosing antiperiodic instead of periodic be across 
the torus; the amplitudes of this relation themselves de- 
pend, in contrast to the 2D case, on the model under 
consideration. 

Note, however, that again in contrast to the 2D case, 
where the influence of boundary conditions on the op- 
erator content has been extensively explored p5| . it is 
theoretically not understood up to now, why using an- 
tiperiodic be in 3D should restore the 2D situation. 



TABLE I. FSS amplitudes of the correlation lengths of the 
Ising, XY, and Heisenberg models on the T 2 x K geometry. 



model 




pbc apbc 


Ising 


Au 
A t 


0.8183(32) 0.23694(80) 
0.2232(16) 0.08661(31) 


A a /A, 

■£e / ■£> <t 


3.666(30) 2.736(13) 
2.7326(16) 


XY 


Aa 

A t 

A a /A t 

/ X (j 


0.75409(59) 0.24113(57) 
0.1899(15) 0.0823(13) 
3.971(32) 2.930(47) 
2.923(7) 


Heisenberg 


A a 
A, 

A a /A e 

Xe 1 


0.72068(34) 0.24462(51) 
0.16966(36) 0.0793(20) 
4.2478(92) 3.085(78) 
3.091(8) 



In view of the total lack of exact results for non-trivial 
3D systems, it seems to us a rewarding challenge for the 
field theorists to explain these results. 
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